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. We report the results of a long observation of the Anomalous X-ray Pulsar 

^ : IRXS J170849-400910 obtained with the BeppoSAX satellite in August 2001. 

<N : The best fit ph ase-averaged spectrum was an absorbed power law plus black- 

body model, with photon index F ~ 2.4 and a black body temperature of 
^ , kThh ~ 0.4 keV. We confirm the presence of significant spectral variations with 

0^ ■ the rotational phase of the pulsar. In the spectrum corresponding to the rising 

part of the pulse we found an absorption-like feature at ~ 8.1 keV (a significance 
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O , of 4cr), most likely due to cyclotron resonant scattering. The centroid energy 

Q ! converts to a magnetic field of 9 x 10^^ G and 1.6 x 10^^ G in the case of electrons 

^ ! and protons, respectively. If confirmed, this would be the first detection of a 

^! cyclotron feature in the spectrum of an anomalous X-ray pulsar. 
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1. INTRODUCTION 

Anomalous X-ray pulsars (AXPs) are characterized by spin periods in the range of 5-12 s, 
steady spin down (~ lO^^"*^ ss^"*^), steep and soft X-ray spectra with luminosities exceding 
by several orders of magnitude their spin-down luminosities (Mereghetti & Stella 1995; van 
Paradijs, Taam & van den Heuvel 1995). All five confirmed AXPs he in the "galactic" plane 
and two (or three), are associated with supernova remnants. AXPs show no evidence for a 
companion and are thus believed to be isolated neutron stars either having extremely strong 
magnetic dipole fields (~ 10^"* — lO-*"^ G; "magnetars"; Duncan & Thompson 1992, Thompson 
& Duncan 1995) or accreting from a residual disk (Li 1999, Alpar 2001, Perna et al. 2000). 
Alternatively, they can accrete from a light companion. IE 1048-5937 and IE 2259+586, 
recently displayed short and intense X-ray bursts (Graviil, Kaspi, & Woods 2002; Kaspi & 
Graviil 2002), strengthening a possible connection between AXPs and Soft 7- Ray Repeaters 
(SGRs; Zhang 2002). For recent rewiews see Israel, Mereghetti & Stella 2002, Mereghetti et 
al. 2002, and references therein. 

IRXS J170849-400910 was discovered with ROSAT (Voges et al. 1996); ~lls pulsa- 
tions were found in its X-ray flux with ASCA (Sugizaki et al. 1997). Early measurements 
suggested that is a fairly stable rotator (Israel et al. 1999). In October 1999 a sudden 
spin-up event occurred, which was interpreted as a glitch (Kaspi et al. 2000). Searches for 
an optical counterpart ruled out the presence of a massive companion (Israel et al. 1999); 
an IR counterpart has been recently proposed (Israel et al. 2003). The diffuse (~ 8') radio 
emission at 1.4 GHz (Gaensler et al. 2001) is due to the supernova remnant G346. 5-0.1, the 
association of which to the AXP is still under debate. There is no evidence of pulsed radio 
emission from the AXP with an upper limit of 70/xJy on the pulsation amplitude (Israel et 
al. 2002). 

Here we report the results of a long BeppoSAX observation of the AXP that confirms 
significant spectral variations with pulse phase (Israel et al. 2001) and shows the presence 
of an absorption-like feature at ~8.1 keV, probably due to cyclotron resonant scattering. 

2. OBSERVATION 

The source was observed by BeppoSAX on 2001 August 17-22 with the imaging Narrow 
Field Instruments (NFI): the Low Energy Concentrator Spectrometer (LEGS: 0.1-lOkeV; 
Parmar et al. 1997; 78.0 ks of effective exposure time) and the Medium Energy Concentrator 
Spectrometer (MEGS: 1.3-1 IkeV; Boella et al. 1997; 200.0 ks exposure). We extracted the 
source photons in the MEGS and LEGS from a circular region of 6 arcmin radius around the 
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position of the AXP. Photons extracted from a sourceless circular region of the same size 
were used for background subtraction. MECS and LEGS photon arrival times were corrected 
to the barycenter of the Solar System. 

3. RESULTS 

3.1. Timing and Spectral Analysis 

In order to obtain a precise estimate of the pulse period, we divided the MECS obser- 
vation in nine time intervals and calculated the pulsation phase for each of them. Fitting 
these phases with a linear function gave a best period of 11.000563 ± 0.000005 s. The folded 
light curve shows an energy-dependent profile (Fig. 1). Specifically, the pulse minimum shifts 
from a phase of ~0.0 in the lowest energy fight curve (0.1-2keV) to ~0.3 in the 6-lOkeV 
fight curve. Correspondingly, the pulsed fraction decreases from ~30% to ~17%. 

We restricted the analysis of the LECS and MECS spectra to the 0.4-4 keV and 1.65- 
10.8 keV energy range, respectively. The spectra were binned so as to have about two 
bins per spectral resolution element (FWHM). Futhermore, the data at the extremes of the 
spectral range given above were further rebinned so as to have at least 20 source events per 
bin (such that minimum chisquare techniques could be reliably used in spectral fitting). 

The spectra were well fit with an absorbed blackbody plus a power law model (see 
Table 1). The best fit of the phase-average spectrum gave a reduced of 0-95 for 298 
degree of freedom (dof) for the following parameters: column density of Nh — (1.36 ± 
0.06) X lO^^cm"^, a blackbody temperature of kTbb = 0.44 ± 0.01 keV (blackbody radius of 
Rbb — 6.6 ± 0.4 km, assuming a distance of 5 kpc) and a photon index of F = 2.40 ± 0.06 (all 
error bars in the text arc 90% confidence). The unabsorbed flux in the 0.5-10 keV range 
was 1.87 X 10^^° ergcm^^ s^^ corresponding to a luminosity of 5.6x10^^ erg s^^ (for a 5 kpc 
distance). In the 0.5-10 keV band the blackbody component accounts for ~30% of the 
total unabsorbed fiux. We also tried other combinations of spectral models (cutoff power 
law plus blackbody or two blackbodies with different temperatures) but they all produced 
larger reduced values. 

3.2. Pulse Phase Spectroscopy 

We adopt for zero phase the minimum of the pulse profile in the 0.2-3 keV folded light 
curve (see Fig. 1). In order to carry out pulse-phase spectroscopy we accumulated spectra 
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in six different pliase intervals. Tfie boundaries of these (0.0, 0.26, 0.4, 0.58, 0.7, 0.84, 1.0) 
were designed so as to sample separately the minimum, rising, maximum, and decaying part 
of the pulse profile, while maintaining a sufficiently good statistics for a detailed spectral 
study (results in Tab. 1, Fig. 2, 3 and 4). A significant variation of the spectral parameters 
with pulse phase was clearly seen (especially for F and Rbb, see Fig. 2; see also Israel et al. 
2001). In all intervals but one, an acceptable fit was obtained with the absorbed power-law 
plus black body model (reduced iii the 0.8-1.1 range); a reduced of 1.2 was instead 
obtained in the 0.4-0.58 phase interval. Specifically, the data were systematically below the 
best fit model in the ~ 7.8-8.4 keV range (see Fig. 4b). We tried to fit three different models: 
a Gaussian, an absorption edge and a cyclotron feature. While the inclusion of a Gaussian or 
an absorption edge did not lead to a significant improvement of the fit, the cyclotron model 
(CYCLABS in the XSPEC package, see Mihara et al. 1990 for details) led to a reduced 
chisquare of 1.1 for 86 dof, corresponding to an F-test probability of 3 x 10~^ {3.1a). 

In order to improve the statistics we added in phase the spectra from the 1999 BeppoSAX 
observation of the source (Israel et al. 2001), increasing the total exposure time to ~ 250 ks. 
For this we used the pulse period determination from the timing solution of a 4 year-long 
i?XTE monitoring (Kaspi et al. 2000) ® . The determination of the zero phase was done again 
by using the pulse profile minimum in the 0.3-2 keV energy range (note that no significant 
shape difference was found with respect to the 2001 observation). We estimate that our 
procedure can introduce an uncertainty of up to 0.01 in the phasing of the two observations 
which is negligible for the aims of our study. 

For each of the six phase intervals, a spectrum was summed together with the cor- 
risponding spectrum from the 2001 observation. The spectral analysis was then repeated. 
In the MEGS and LEGS spectra from the 0.4-0.58 phase interval the inclusion of a Gaus- 
sian (£'g = 8.3keV, cr = 0.4keV and x^ = l-l) and an absorption edge (£'e = 7.7keV, t = 0.4 
and x^ = 0.9) to fit the feature at ~ 8.1 keV resulted in an improvement of the chisquare, 
which converted to a single-trial F-test probability of 0.18 (~ la) and 4 x 10~^ (~ 3.7a), 
respectively. A much more significant improvement was obtained by adding instead a Res- 
onant Cyclotron Feature (RCF) model; the F-test probability in this case was 1.8 x 10^^ 
corresponding to a single trial significance of 4.5(7 or Aa after correction for the six spectra 
that we analysed (see Fig. 4 and Tab 1). 



^Note that we could not use for the 2001 dataset the post-"gUtch" timing (Gavriil & Kaspi 2002) since 
the period extrapolated from their parameters significantly differs from that derived here (AP ~ 1.5 x 10~^ 
s) 
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4. DISCUSSION 

During a BeppoSAX study of IRXS J170849— 400910 wc discovered an absorption-like 
feature at an energy of ~8.1keV in a pulse phase interval corresponding to the rising part 
of the ~ lis pulse. This feature was best fit by a resonant cyclotron feature model with a 
centroid energy of ~8.1keV and an equivalent width of ~460eV. 

The detection of an RCF in a specific pulse-phase interval and superposed to an X-ray 
continuum that varies with the pulse phase is reminiscent of the behaviour seen in standard 
accreting pulsars in X-ray binaries (see Wheaton et al. 1979, Santangelo et al. 1999). If 
interpreted as an electron resonant feature at the base of the accretion column, the feature 
at ~8.1 keV imphes a neutron star surface magnetic field of ~ 9.2 x 10^^ Gauss (using 
a gravitational redshift z=0.3) . This value is just sfightly lower than that measured for 
electron RCFs in typical accreting X-ray pulsars (see Fig. 5); more interestingly it is close 
to that required by models for AXPs which involve residual disk accretion in the spin-down 
regime (Mereghetti & Stella 1995; Alpar 2001; Perna et al. 2001 ). In this context, one can 
solve the torque equation (see e.g. Eq. 11.35 in Henrichs 1983) by exploiting the measured 
value of P and range of accretion luminosity derived from plausible distances (5-10 kpc). 
The surface magnetic (dipole) field obtained in this way is 0.6-1.1x10^^ G (corrisponding to 
a fastness parameter range of cUg — 0.57 — 0.54, a typical value for the spin-down accretion 
regime; see Ghosh & Lamb 1979 and Henries 1983). 

The agreement of this estimate with the magnetic field inferred from the electron RCF 
interpretation is intriguing, especially in consideration of the other analogies with the pulse- 
phase spectral dependence of conventional accreting X-ray pulsars. By contrast, if an electron 
RCF arose somehow at the polar caps of a rotation powered pulsar, a B-field strength of 
9.2x10^^ G would be in the range of many radio pulsars and yet much lower than that 
required to spin-down at the observed rate through magnetic dipole radiation (~ 5 x 10^^ G, 
indeed this was one of the motivations for magnetar model, see below) . 

Alternatively the RCF might be due to protons. For the magnetic field strengths forseen 
in the "magnetar" scenario, proton cyclotron features (if any) are expected to lie in the 
classical X-ray band (0.1-lOkeV; Zane et al. 2001; Lai & Ho 2002). A proton RCF feature 
at ~ 8.1keV would correspond to surface field of 1.6x10^^ G (z=0.3). The fact that this 
value is ~ 3 times higher than the surface field derived from the usual magnetic dipole spin 
down formula should not be of concern. According to the magnetar model, the magnetic field 
at the star surface and its vicinity is dominated by higher order multipole field components 
(Thompson & Duncan 1995). At large radii the dipole component, responsible for the secular 
spin-down, dominates. It is thus expected that a proton RCF feature, sampling the (total) 
surface magnetic field strength, provides a higher value than the mere dipole component. 
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There is a clear correlation between the width and centroid energy of the electron RCFs 
in accreting X-ray pulsars (see Fig. 5 extending the results of Orlandini & Dal Fiume 2001). 
The values from IRXS J170849-400910 and SGR 1806-20 are in good agreement with such 
a relation. The modest range of width to centroid energy ratio implied by this indicates 
that magnetic field geometry effects at the neutron star surface likely dominate the RCF 
width (on the contrary temperature and particle mass would alter this ratio). This, in turn, 
suggests that similar (relative) ranges of surface magnetic field strength are "sampled" by 
RCFs in accreting X-ray pulsar and, by extension, RCFs in AXPs and SGRs. 

Other interpretations of the feature at ~ 8.1 kcV appear less likely. Firstly, fitting an 
edge or line due to photo-electric absorption provides a less pronounced improvement of the 
fit than the RCF model. Secondly, an edge by iron at a sufficiently large distance from the 
neutron star that energy shifts are negligible would require a high overabundance of this 
element and intermediate ionisation stages (such as C-like iron) . Yet it has long been known 
that the photoionisation equilibrium of such a plasma is unstable (Krolik & Kallman 1984; 
Nagase 1989). The energy of an ion feature forming in the neutron star atmosphere would be 
drastically altered by magnetic field effects (see Mori & Hailey 2002 and references therein). 
In this and the above interpretations, however, it would also be difficult to explain why an 
ion feature is observed only over a restricted range of pulse phases. Exploring in detail these 
possibihties is beyond the scope of this letter. 

In conclusion, we found an absorption-like feature in the BeppoSAX X-ray spectrum 
of IRXS J170849— 400910 taken during the rising phase of the ~lls pulse, which is best 
fit by a cyclotron resonant scattering model. If this interpretation is correct, the centroid 
energy translates into a magnetic field strength of ~1.6xl0^^ G and ~9.2xl0^^ G depending 
on whether protons or electrons, respectively, are responsible for the feature. 

Nanda Rea acknowledges useful discussions with S. Dall'Osso. This work was partially 
supported through ASI and COFIN 2000 grants. 
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Table 1: Best fit spectral parameters from selected pulse- phase intervals: that showing feature 
at ~8.1kcV (0.4-0.58) and those characterized by the lowest and the highest power law 
photon index T (see also Fig. 3). All fluxes are unabsorbed and calculated in 0.5-10 keV 
range; uncertainties are 90% confidence. 



Parameters 


0.4-0.58 


Phase Intervals 
0.84-1.0 


0.26-0.4 


Nh{x10^^ cm-2) 


1.38 ±0.01 


1.1 ±0.1 


1.6 ±0.3 


r 


2.6 ±0.1 


2.0 ±0.1 


2.9 ±0.1 


PL flux (xlO~^^ergcm~^ s~^) 


21. 


8.3 


24.8 


(keV) 


0.46 ± 0.01 


0.46 ± 0.02 


0.42 ± 0.01 


BB flux (xlO-^^ergcm-^s-^) 


9.1 


5.1 


6.1 


Rbb (d = 5kpc; km) 


6.7 ± 0.4 


5.0 ±0.6 


7.2 ± 0.6 


Ecyclabs (keV) 


8.1 ±0.1 






Line Width ( keV) 


0.2 ±0.1 






Line Depth ( keV) 


0.8 ±0.4 






xVd.o.f. 


1.09 


1.01 


0.98 


Flux (xlO~^^ ergcm"^ s~^) 


30.0 


7.3 


7.2 
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Fig. 1.— MECS light curves of IRXS J170849-400910 folded at the best spin period (two 

pulse cycles are shown) for six energy bands: (a) 0.1-2 keV; (b) 2-3 keV; (c) 3-4 keV; (d) 

1 1 
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Fig. 2. — Spectral parameter variability from pulse-phase spectroscopy using an absorbed 
black body plus power-law model. Filled squares represent the spectral parameters after the 
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Fig. 3. — LEGS and MEGS spectra in the two phase intervals characterised by the highest 
and lowest value of the power law photon index F; see also Table 1. 
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Fig. 4. — MECS and LEGS spectra from the 0.4-0.58 phase interval fit with the "standard 
model" (the sum of a blackbody and power law with absorption) plus a cyclotron line. 
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Fig. 5. — Line width vs. centroid energy from a sample of accreting X-ray pulsars with 
electron RCFs (Orlandini & Dal Fiume 2001), IRXS J170849-400910 (this Letter) and 
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